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Abstract — Max-product belief propagation is a local, iterative 
algorithm to find the mode/MAP estimate of a probability distri- 
bution. While it has been successfully employed in a wide variety 
of applications, there are relatively few theoretical guarantees 
of convergence and correctness for general loopy graphs that 
may have many short cycles. Of these, even fewer provide exact 
"necessary and sufficient" characterizations. 

In this paper we investigate the problem of using max-product 
to find the maximum weight matching in an arbitrary graph with 
edge weights. This is done by first constructing a probability 
distribution whose mode corresponds to the optimal matching, 
and then running max-product. Weighted matching can also be 
posed as an integer program, for which there is an LP relaxation. 
This relaxation is not always tight. In this paper we show that 

1) If the LP relaxation is tight, then max-product always 
converges, and that too to the correct answer. 

2) If the LP relaxation is loose, then max-product does not 
converge. 

This provides an exact, data-dependent characterization of max- 
product performance, and a precise connection to LP relaxation, 
which is a well-studied optimization technique. Also, since LP 
relaxation is known to be tight for bipartite graphs, our results 
generalize other recent results on using max-product to find 
weighted matchings in bipartite graphs. 

I. Introduction 

Message-passing algorithms, like Belief Propagation and 
its variants and generalizations, have been shown empir- 
ically to be very effective in solving many instances of 
hard/computationally intensive problems in a wide range of 
fields. These algorithms were originally designed for exact 
inference (i.e. calculation of marginals/max-marginals) in tree- 
structured probability distributions. Their application to gen- 
eral graphs involves replicating their iterative local update 
rules on the general graph. In this case however, there are 
no guarantees of either convergence or correctness in general. 

Understanding and characterizing the performance of 
message-passing algorithms in general graphs remains an 
active research area. [1,2] show correctness for graphs with 
at most one cycle. [3,4] show that for gaussian problems 
the sum-product algorithm finds the correct means upon con- 
vergence, but does not always find the correct variances. [5, 
6] show asymptotic correctness for random graphs associated 
with decoding. [7] shows that if max-product converges, then 
it is optimal in a relatively large "local" neighborhood. 

In this paper we consider the problem of using max-product 
to find the maximum weight matching in an arbitrary graph 
with arbitrary edge weights. This problem can be formulated 



as an integer program, which has a natural LP relaxation. In 
this paper we prove the following 

1) If the LP relaxation is tight, then max-product always 
converges, and that too to the correct answer. 

2) If the LP relaxation is loose, then max-product does not 
converge. 

Bayati, Shah and Sharma [8] were the first to investigate 
max-product for the weighted matching problem. They showed 
that if the graph is bipartite then max-product always con- 
verges to the correct answer Recently, this result has been 
extended to 6-matchings on bipartite graphs [9]. Since the LP 
relaxation is always tight for bipartite graphs, the first part 
of our results recover their results and can be viewed as the 
correct generalization to arbitrary graphs, since in this case 
the tightness is a function of structure as well as weights. 

We would like to point out three features of our work: 

1) It provides a necessary and sufficient condition for con- 
vergnce of max-product in arbitrary problem instances. 
There are very few non-trivial classes of problems for 
which there is such a tight characterization of message- 
passing performance. 

2) The characterization is data dependent: it is decided 
based not only on the graph structure but also on the 
weights of the particular instance. 

3) Tightness of LP relaxations is well-studied for broad 
classes of problems, making this chracterization promis- 
ing in terms of both understanding and development of 
new algorithms. 

Relations, similarities and comparisons between max-product 
and linear programming have been used/mentioned by several 
authors [10-12], and an exact characterization of this relation- 
ship in general remains an interesting endeavor. In particular, 
it would be interesting to investigate the implications of these 
results as regards elucidating the relationship between iterative 
decoding of channel codes and LP decoding [13]. 

II. Weighted Matching and its LP Relaxation 

A matching in a graph is a set of edges such that no two 
edges in the set are incident on the same node. Given a graph 
G = {V,E), with non-negative weights We on the edges e S 
E, the weighted matching problem is to find the matching M* 
whose edges have the highest total weight. In this paper we 
find it convenient to refer to edges both as e G -E and as 
where i,j e V. 



Weighted matching can be written as the following integer 
program (IP): 

max WeXe 

s.t. ^ < 1 for alH e (1) 

Xe G {0, 1} for all e e 

The LP relaxation of the above problem is to replace the 
constraint Xe £ {0,1} with the constraint Xe > 0. This 
relaxation is in general not tight, i.e. there might exist non- 
integer solutions with strictly higher value than any integral 
solution. It is known however that the LP relaxation is always 
tight for bipartite graphs: no matter what the edge weights, the 
bipartite-ness ensures tightness of the LP relaxation. If a graph 
is not bipartite, the tightness of the LP relaxation will depend 
on the edge weights: the same graph may have tightness for 
one set of weights and looseness for another set. 

The dual of the above linear program is the vertex cover 
problem: minimize the total of the weights Zi that need to be 
placed on nodes so as to "cover" the edge weights: (DP) 

min Zi 

s.t. Wij < Zi + Zj for all (i, j) e E 

Zi > for all i 

Lemma 1 (complimentary slackness): When the LP relax- 
ation is tight, the optimal matching M* and the optimal dual 
variables z and satisfy the following properties: 

1) if («, j) G M* then Wij — Zi + zj 

2) if («, j) ^ M* then Wij < Zi + Zj 

3) if no edge in M* is incident on node i, then Zi = Q 

4) Zi < maxe We for all i 

III. Background on the Max-Product Algorithm 

The factor graph [14] of a probability distribution represents 
the conditional independencies of the distribution. The Max- 
Product (MP) algorithm is a simple, local, iterative message 
passing algorithm that can be used (in an attempt) to find 
the mode/MAP estimate of a probability distribution. Nodes 
and factors pass messages to each other, and nodes maintain 
"beliefs", which represent the max-marginals. When max- 
product is applied to problems involving general "loopy" 
graphs, one of the following three scenarios may result: 

1) The algorithm may not converge. 

2) The algorithm may converge, but to an incorrect answer 

3) The algorithm may converge to the correct answer 

As has been mentioned, here has been siginifcant work at- 
tempting to understand the properties of MP for loopy graphs. 
For the results in this paper, we will use the following two 
insights: 

1) At any time, the belief of the max-product algorithm 
for a given variable corresponds to the belief at the root 
of the corresponding computation tree distribution [2] 
associated with that variable at that time. We describe 



what this computation tree distribution corresponds to 
for the weighted matching problem in the next section. 
2) If max-product does converge, the resulting beliefs are 
optimal in a large "local" neighborhood [7]: let x be 
the assignment as given by the converged max-product 
and X be any other assignment. If the variables assigned 
different values in x and x form an induced graph 
containing at most one cycle in each component, then 
p{x) > p{x). 

IV. Max-Product for Weighted Matching 

The problem of finding M* can be formulated as the prob- 
lem of finding the mode of a suitably (artifically) constructed 
probability distribution p. In fact, there are in general several 
ways to construct this distribution for the same instance of a 
graph G. We now present one constructiorQ. 

Associate a binary variable Xe £ {0, 1} with each edge 
e e E, and let 

Pi-) = ini{i:....,-.<i} n^'"^^^ (2) 

iev eeE 
Here J\f{i) represents the neighborhood of node i in G, and Z 
is a normalizing constant. The variable Xe can be interpreted 
as follows: = 1 indicates that e G M*, while Xe = 
indicates e ^ M* . The term > enforces the 

cosntraint that of the edges incident to node i, at most one 
can be assigned the value "1". Thus, it is easy to see that 
p{x) > if and only if the edges with x^ — I constitute a 
matching in G. Furthermore, the mode of p corresponds to the 
max- weight matching M*. 

The factor graph max-product involves messages between 
variables and factors. In our case the variables are the edges 
G E, and the factors are nodes i E V. Thus at any 
time t there will be messages m^_^^_^ from node (factor) i 
to edge (variable) as well as messages m*. ^j^^^. Each 

message will be a length-two vector of real numbers, indexed 
by and 1 . The message update rules can be simplified to the 
following: 

k(iN(i)-j 
= max{ [] m\^,^)^M . 

keJ\f{i)-j 

Also, at every time each edge (variable) maintains a belief 
vector 6*j as follows: 

'This construction is different from the one in [8], which had a pairwise 
model with variables corresponding to nodes in the graph. However, the results 
of this paper continue to hold when the construciton in [8] is modified to be 
applicable to general graphs 



The p defined above can be used to find M* as follows: first 
run max-product. At any time t and for each edge e there will 
be two beliefs 6* [0] and bl[l]. If max-product converges, assign 
to each variable the value (i.e. "0" or "1") that corresponds 
to the stronger belief. Then, declare the set of all edges set to 
"1" to be the max-product output. 

A. The Computation Tree for Weighted Matching 

Our proofs rely on the computation tree interpretation [2, 
15] of the Max-product beliefs. We now describe this inter- 
pretation when max-product is applied to p as given in (|2|. 

For an edge e let Te{k) be the full depth-k computation tree 
rooted at e. This is generated recursively: take Te{k — 1) and 
to each leaf v add as children a copy of each of the neighbors 
of V in G, except for the unique neighbor of v which is already 
present in Te{k — 1). Also, each new edge has the same weight 
as its copy in the original G. The recursion is started with 
the single-edge tree Te(l) = e, both of whose endpoints are 
leaves. This initial edge is the root of Tg. 

Consider now the "full synchronous" max-product, where at 
each time every message in the network is updated. In this case 
the computation tree Te{k) for edge e at time k will be Te{k). 
Alternatively, max-product may be executed asynchronously 
with only a subset of the messages updated in every time slot. 
In this case Te{k) will be a sub-tree of Te{k). In either case, 
the computation tree interpretation states at time k we have 
&g[l] > 6g[0] if and only if the root of Tf,{k) is a member of 
a max-weight matching on the tree Te{k). 

The figure below shows an example where on the left is 
G: the four-cycle abed and the chord ac, with a matching 
M — {(a, 6), (c, d)} depicted in bold. On the right is the 
computation tree Ti^a.b){'^) which is the full tree of depth 4 
rooted at edge (a, b). The bold edges depict the projection Mt 
of M onto j')(4): an edge e in the tree is in Mt if and 
only if its copy in G is in M. 
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Lemma 2: Let M be a matching in G and Te(fc) be a 
computation tree. Let Mt be the set of all copies in Te{k) 
of all edges in M. Then, Mt is a matching in Te{k). Also, if 
M is maximal in G, Mt is maximal in Tg. 

Of course Tg will also contain other matchings that are not 
projections of matchings in G. Finally, we say that a (possibly 
not full) tree Te(fc) Is/mW upto depth ki if the full tree Te{ki) 
is contained in Te{k). 



V. Equivalence of Max-Product and LP 
Relaxation 

We are now ready to prove the main result of this paper: 
the equivalence of Max-Product and LP Relaxation. Before 
we proceed, we define the following terms 

1) We say that the LP relaxation is tight if the linear 
program (LP) obtained by relaxing the integer program 
([TJ has a unique optimal solution at which all values 
are either or 1. 

2) We say that max-product converges by step k if the 
variable assignments (0 or 1) that maximize the beliefs 
at each node remain constant once the associated com- 
putation tree is full up to depth at least k. Note that this 
includes both synchronous and asynchronous message 
updates. We say that max-product converges if there 
exists some k < oo such that max-product converges 
by step k. Finally, we say that max product converges 
to the correct answer if the beliefs fee at convergence 
are such that fee[l] > fee[0] if and only if e G M*, and 
fee[l] < fee[0] if and only if e ^ M* 

We also need to make some uniqueness assumptions. It is 
well-recognized that max-product may perform poorly in the 
presence of multiple optima, and that characterizing perfor- 
mance in this case is hard. For the rest of this paper we will 
assume the following: 
Al M* is the unique optimal matching. 
A2 The linear program always has a unique optimal solu- 
tion. Note that this can be fractional, but it has to be 
unique. 

A. Max-product is as Powerful as LP Relaxation 

In this section we prove that if the LP relaxation is tight then 
Max-Product converges to the correct answer Recall that when 
the LP is tight, part 2 of Lemma [T] says that if ^ M* 
then Wij < Zi + Zj. The uniqueness assumptions Al-2 further 
imply that the inequality is strict: Wij < Zi + zj. Another way 
of saying this is that there exists an e > such that 

Wij < Zi + Zj — e for all ^ M* (3) 

Theorem 1: Consider a weighted graph G for which the LP 
relaxation is tight. Then max-product converges to the correct 
answer by step where Wmax = maxe is the weight 

of the heaviest edge, and e satisfies (|3]l. 
Proof: 

Let M* be the optimal matching on G. For max-product to 
be convergent and correct, we need that fe*[l] > fe*[0] for all 
e e M* and fe*[l] < fe*[0] for all e ^ M*, and for all t such 
that Te(t) is full upto depth ^liim- 

So suppose that for such a t there exists an e ^ M* such 
that fe*[l] > fe*[0]. Then, there exists a matching M in Te{t) 
such that (a) the root e e M, and (b) M has the largest weight 
among matchings on Te{t). Let AI^ be the set of all edges 
in Te{t) that are copies of edges in M*. By lemma |2] M^ 
is a maximal matching on Te{t). Also, the root e ^ AI^ by 
assumption. 



The symmetric difference M^AM consists of disjoint 
alternating paths in Te{t): each path will have every alternate 
edge in and all other edges in M. Let P be the path 
that contains the root e. We now show that w{P n M^) > 

wiPnM). 

Recall that the optimal dual solution assigns to each node i 
in G a "dual value" Zi > 0. Associate now with each node in 
Te (t) the dual value of its copy in G. Then, by Lemma [T] we 
have that Wij = Zi + Zj for each (i, j) G P n M^. Suppose 
now that neither endpoint of P is a leaf of Te(i). In this case, 
we have 

ieP 

On the other hand, we know that (|3]l holds for each edge in 
P n M. Adding these up gives 

u)(P nAf) < ^z, -e|PnAf| 

By assumption, the root e G P fl M, so |P n Af | > 1 and 
hence w(P n M^) > w{P D M) when no endpoints of P are 
leaves. 

Suppose now that exactly one of the endpoints w of P is a 
leaf of Te{t). In this case, we have that 

w{P n M^) > 

ieP iGP 

where the last inequality follows from part 4 of Lemma[T] Also, 
Te{t) is assumed to be full up to depth k, so this implies that 
P n Af| > |. This means that 

«;(PnA/) < ^z, -e^ 

ieP 

Now, since k > ^"''J'"" , this impHes that w(PnA/^) > w{Pn 
Af). The final case, where both endpoints of P are leaves, 
works out in the same way, except that now |Pn A/| > k and 

W(P n M}) > J^teP - '^Wmax- 

Thus, in any case, we have that w{P n M^) > w{P fl Af). 
Consider now the set of edges Af - (P n Af ) + (P n Af^). 
This set forms a matching on Te{t), and has higher weight 
than M. This contradicts the choice of Af , and so establishes 
that < 6*[0] for all e ^ M*. A similar contradiction 

argument can be used to establish that > 6*[0] for all 
e G Af*. This completes the proof. ■ 

B. LP Relaxation is as Powerful as Max-product 

In this section we prove that if the LP relaxation is loose 
then max-product does not converge to the correct answer. 
Before we do so however, we note that this implies a stronger 
result: that when LP is loose then in fact max-product does 
not converge at all. 

Lemma 3: Consider the distribution p{x) as given in (|2|i. If 
Max-Product converges, then its output exactly corresponds to 
the true optimal matching Af*. 

The proof of this lemma uses the "local optimality" result 
of Weiss and Freeman [7]. In particular, for p it turns out that 
local optimality implies global optimality. This means that it 



is not possible for max-product to converge to an incorrect 
answer: it will either not converge at all, or converge to Af *. 
We do not use this explicitly in the proofs below, but it 
strengthens the results as mentioned above. 

We now proceed with showing that max-product does not 
converge to the correct Af* when LP is loose. As a first 
step, we need a combinatorial characterization of when the 
LP relaxation is loose. We now make some definitions. We 
say that a node v is saturated by a matching Af if there exists 
an edge e G Af that is incident to v. 

A blossom with respect to a matching Af is an odd cycle 
C with ed ges m AfS Note that a blossom has a unique 

base: a node not saturated by any edge in C n Af . A stemmed 
blossom Bi (w.r.t Af ) is a blossom C, along with an alternating 
path (stem) P that starts at the base of C, and starts with an 
edge in Af . Also, P should be such that the set Af — (P n 
Af ) + (P — Af ) remains a matching in G. 

A bad stemmed blossom is one in which the edge weights 
satisfy 

w(C n Af ) + 2w(P n Af ) < wiC - M) + 2w{P - M) 

Note that it may well be the case that |P| = 0, in which case 
Bi is just an odd cycle. The following is an example of a 
bad stemmed blossom. The bold edges are the ones in Af , the 
numbers denote the weights of the corresponding edges, and 
the last node i has no edge of Af incident on it. The blossom 
C in this case is the cycle abode, and node c is its base. The 
path/stem P is cfghi. 



b 




A blossom pair B2 is two blossoms Ci and C2 and an 
alternating path P between the bases of the two blossoms 
such that P begins and ends with edges in Af . A bad blossom 
pair is one in which the edge weights satisfy 

w(Ci n Af ) + w(C2 n Af ) + 2w{P n Af ) 

< w{Gi - Af ) + w{C2 - Af ) + 2w(P - Af ) 
The following is an example of a bad blossom pair. 




-Blossoms were first defined in [16], which also provided the first efficient 
algorithm for weighted matching in arbitrary graphs. 



The following proposition provides a combinatorial charac- 
terization of when the LP relaxation is loose, and is crucial to 
the proof of the subsequent theorem. 

Proposition 1: If the LP relaxation is loose, then there 
exists a bad stemmed blossom, or a bad blossom pair, with 
respect to the optimal matching M* . 

Proof: In appendix. 

We use the presence of these "bad" subgraphs in G to show 
that max-product does not converge to the correct answer. 
Before we do so, we need one additional lemma. This states 
that if max-product converges by step k to some matching M 
on G, then the optimal matching Mt on the computation tree 
looks like M in the neighborhood of the root. 

Lemma 4: Suppose max-product converges to a matching 
M in G by step k. Consider any edge e, some m > 1 and a 
corresponding computation tree Tg which is full up to depth 
k + m. Let Mt be the max- weight matching on the tree. Then, 
for any edge f ^ that is within distance m of the root e, 
/ G Mt if and only if its copy /i in G is such that /i G M. 

Note that the above lemma also applies to the root e of the 
tree. We are now ready to state and prove the main result of 
this section. Recall that the belief be on an edge at convergence 
is incorrect if either e £ M* but be[0] > be[l], or e ^ M* but 

be[l] > be[0]. 

Theorem 2: Consider a weighted graph G for which the 
LP relaxation is loose. Then, the max-product beliefs do not 
converge to the correct M*: for any given k, there exists a 
ki > k and computation trees Te,e £ E such that each Tg is 
full upto depth fci, but the beliefs on some of the edges are 
incorrect. Lemma |3] further implies that in fact in this case 
max-product does not converge at all. 
Proof: 

Let M* be the max-weight matching on G. Since the 
LP relaxation is loose, by Prop. [T] there exists either a bad 
stemmed blossom or a bad blossom pair w.r.t. M* . Suppose 
first that it contains a bad stemmed blossom Bi, and consider 
some e G C n M* that is in the "blossom" part of Bi (as 
opposed to the stem) and also in M* . From the two nodes of 
e, make maximal alternating paths Pi and P2 that remain in 
Bi and start out in opposite directions on C. For the stemmed 
blossom example above, if e is the edge (a, b) then the two 
paths will be bcfghi and aedcfghi. 

Let di = w{Pi - M*) - w{Pi n M*), and similarly ^2 
for P2 ■ di represents the change in the weight of the matching 
if each edge in Pi were "switched", i.e. their membership in 
the matching was reversed from its original value. It is easy 
to see that 

di+d2-w{e) = w{C - M*) + 2w{P - M*) 

~w{C n M*) - 2w{P n M*) 

By assumption Bi is a bad blossom and hence we have that 
di+d2- w{e) > 0. 

Suppose max-product converges to M* by step k. Consider 
now the computation tree Te which is full upto depth k + 
\V\, where is the number of nodes in G. Let Mt be the 



max-weight matching on Te- Lemma |4] implies that Mt will 
be a projection of M* in a distance-|y| neighborhood of the 
root. Also, starting from the root e, each of Pi and P2 will 
have a unique copy, say Ri and R2 respectively, in T^, with 
|i?i|, |i?2| < \V\. Since Pi and P2 are alternating w.r.t. M*, 
it follows that Ri and R2 will be alternating with respect to 
Mt- Also, the set 5 = i?i U e U i?2 forms an alternating path 
on Te with respect to Mt, and this begins and ends in nodes 
unsaturated by Mt- Thus, AIt can be augmented by this path: 
the set Mt ~ {S n Mt) + {S ~ Mt) will be a matching on 

Te- 

Also, the weight gain from doing this augmentation will be 
exactly di + d2 — ^(e), which we know is strictly positive. 
Thus, this shows that Mt is not the optimal matching on 
Te, which contradicts the choice of Mt- This means that 
our assumption about max-product convergence to M* is 
incorrect. 

Thus, we see that if there exists a bad stemmed blossom 
w.rt. M* in G then max-product does not converge to M*. 
A similar argument holds for the case of a bad blossom pair 
B2, except that instead of paths Pi and P2 above we now 
have to look at alternating walks Wi and W2 that live in B2 
and are long enough. These walks can then be mapped to an 
augmenting path on Tg which strictly improves Mt, leading 
to a contradiction as was seen in the case of the paths Pi and 
P2- This completes the proof. ■ 

VI. Discussion 

The results of this paper can be generalized to the case 
of perfect matchings, 6-matchings and perfect 6-matchings in 
general graphs, where similar results hold. In this paper max- 
product is shown to be as powerful as LP relaxation, but it 
would be more interesting to outline a direct operational link 
between max-product and a linear programming algorithm. As 
an example, [8] shows that for bipartite matching max-product 
has an operational correspondance with the auction algorithm 
[17]. Also, the form of the message update equations suggests 
that it can be implemented via an equivalent message passing 
update rule between just the nodes of the graph G, instead of 
having messages go from nodes to edges and vice versa. 

More generally, it would be interesting to see if the ideas 
presented in this paper could be used/genealized to show con- 
nections between linear programming and belief propagation 
in other applications. 
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Appendix 

Proof of Proposition [T] 

We now show that if the LP relaxation is loose then there 
exists in the graph either a bad stemmed blossom or a bad 



blossom pair, with respect to the optimal matching M*. Let 
X be the optimal (fractional) solution to the LP relaxation. 

Let E' be the set of all edges e such that either (a) e G M*, 
or (b) e ^ M* and > 0. Then, E' will contain at least one 
edge e ^ M*, because if all e ^ M* had Xe = then the 
LP would be tight. Let G' = {V,E') be the subgraph of G 
having only the edges in E'. An cycle augmentation is any 
even cycle in which every alternate edge is in M*. A path 
augmentation is any path in which every alternate edge is in 
M*, and which begins and ends in nodes unsaturated by M*. 
For any augmentation A, we have that M* -{AnM*) + {A- 
M*) is also a matching in G'. Thus, if M* is the unique max- 
weight matching it has to be that w{A fi M*) > w{A — M*). 

Lemma 5: G' cannot contain any augmentations: cycles or 
paths. 

Proof: Let A be an augmentation in G'. By assumption, 
a;e > for all e € ^ — M*, which imphes that Xg < 1 for 
all e G ^ n M*. Thus, there exists some e > such that 
decreasing each Xe,e £ A ~ M* by e and increasing each 
Xe,e € An M* by e represents a vahd new feasible point 
for the LP. The weight of this new point exceeds the weight 
of X by e{xv{A D M*) - ■u}{A - M*)) > 0. However this 
contradicts the optimality of x, and thus G' cannot contain 
any augmentation. ■ 

Let S be the longest alternating sequence of edges in G", 
and let vi and V2 be its endpoints. By the lemma above, both 
cannot be unsaturated. We say that vi or V2 is a saturated leaf 
if it is saturated by M* and there exist no edges in G' — M* 
incident on it. Also, note that an endpoint is saturated if and 
only if its corresponding edge in S is also in M*. 

The fact that S is the longest sequence means that it cannot 
be extended further beyond vi and V2. This implies that one 
of the following cases must occur: 

1) Both vi and V2 are both saturated leaves 

In this case, the constraints at vi and are loose. So, 
there exists an e such that if all aig, e e S — M* are 
decreased by e and all Xe,e G S D M* are increased 
by e then the new solution remains feasible. This new 
solution will have strictly higher weight than x, which 
is a contradiction. Thus this case cannot occur. 

2) vi is a saturated leaf and V2 is unsaturated. 

An e-perturbation argument like the one above can be 
used to show that this case too cannot occur. 

3) vi is saturated by M*. but is not a leaf. V2 is either 
unsaturated, or a saturated leaf. 

Since S cannot be extended, it has to be that all edges in 
G' — M* incident to vi have other endpoints in S. Let e 
be one such edge. Then, en5 forms a stemmed blossom: 
the resulting cycle has to be odd, and the remaining part 
of S will be a stem whose endpoint is V2. Note that in 
this case it has to be that the constraint at V2 is loose. 

4) Both vi and V2 are saturated by M*, but are not leaves. 
Applying the above blossom argument to both vi and 
V2 yields the existence of a blossom pair. 

Thus if the LP relaxation is loose then there exists a 
stemmed blossom or a blossom pair. Now all that is remaining 



to show is that they are "bad". Let Bi be a stemmed blossom 
in G', consisting of blossom C and stem P. Then, there exists 
some e > such that if Xe,e G C H M* is increased by e, 
Xe, e S C — M* is decreased by e, Xe, e G PnM* is increased 
by 2e, and Xe,e € C — M* is decreased by 2e, then the new 
solution remains feasible for the LP. Also, the new solution 
weighs 

e [w{C n M*) + 2w{P n M*) - w{G - M*) - 2w{P - M*)] 

more than x. For x to be the unique optimal of the LP, this has 
to be strictly negative and thus any stemmed blossom Bi is 
bad. A similar argument shows that any blossom pair is bad. 
This finishes the proof of the proposition. ■ 
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